cosmological perturbation theory 

Revisited 



Claes Uggla* 

^ ■ Department of Physics, 

O ; University of Karlstad, S-651 88 Karlstad, Sweden 

q' 

^ : John WAiNWRiGHTt 



X 



Department of Applied Mathematics, 
University of Waterloo, Waterloo, ON, N2L 3G1, Canada 



<N : January 13, 2013 

cn 

o 
in 

(N 

o 



Abstract 



Increasingly accurate observations are driving theoretical cosmology to- 
ward the use of more sophisticated descriptions of matter and the study of 
nonlinear perturbations of Friedmann-Lemaitre cosmologies, whose governing 
equations are notoriously complicated. Our goal in this paper is to formu- 
late the governing equations for linear perturbation theory in a particularly 
■ simple and concise form in order to facilitate the extension to nonlinear per- 

turbations. Our approach has several novel features. We show that the use 
of so-called intrinsic gauge invariants has two advantages. It naturally leads 
to: (i) a physically motivated choice of a gauge invariant associated with the 
matter density, and (ii) two distinct and complementary ways of formulating 
the evolution equations for scalar perturbations, associated with the work of 
Bardeen and of Kodama and Sasaki. In the first case the perturbed Einstein 
tensor gives rise to a second order (in time) linear differential operator, and 
in the second case to a pair of coupled first order (in time) linear differential 
operators. These operators are of fundamental importance in cosmological 
perturbation theory, since they provide the leading order terms in the govern- 
ing equations for nonlinear perturbations. 
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1 Introduction 

The analysis of linear perturbations of Friedmann-Lemaitre (FL) cosmologies was 
initiated by Lifshitz (1946) in a paper of far-reaching importance. Working in the 
so-called synchronous gauge, this paper showed that an arbitrary linear perturbation 
can be written as the sum of three modes, a scalar mode that describes perturbations 
in the matter density, a vector mode that describes vorticity and a tensor mode that 
describes gravitational waves. For many years, however, the theory was plagued by 
gauge problems, i.e. by the fact that the behaviour of the scalar mode depends 
significantly on the choice of gauge. A major step in alleviating this difficulty was 
taken by Bardeen (1980), who reformulated the linearized Einstein field equations 
in terms of a set of gauge-invariant variables, as an alternative to the traditional use 
of the synchronous gauge. Central to Bardeen's paper are two gauge-invariant equa- 
tions that govern the behaviour of scalar perturbations. The first of these governs 
the evolution in time of a gauge-invariant gravitational [i.e. metric) potential and 
the second determines a gauge-invariant perturbation of the matter density in terms 
of the spatial Laplacian of the gravitational potential. Since this potential continues 
to play a central role in the study of scalar perturbations, it seems appropriate to 
refer to it as the Bardeen potential. Bardeen's paper makes clear, however, that 
there is no unique way of constructing gauge-invariant variables. 

From our perspective, one drawback of Bardeen's paper is that he performs a 
harmonic decomposition of the variables ab initio, with the result that the mathe- 
matical structure of the governing equations is somewhat obscured. In a subsequent 
paper, Brandenberger, Khan and Press (1983) address this deficiency by giving a new 
derivation of Bardeen's gauge-invariant equations. They do not perform a harmonic 
decomposition, with the result that their evolution equation is a partial differential 
equation rather than an ordinary differential equation as in Bardeen's paper. How- 
ever, unlike Bardeen they restrict consideration to a spatially flat Robertson- Walker 
(RW) background^ 

In subsequent developments the status of the Bardeen potential was further en- 
hanced by the appearance of the major review paper by Mukhanov et al (1992), 
which contains a simplifled derivation of the Bardeen potential and the evolution 
equation for scalar perturbations, without performing a harmonic decomposition. 
However, the treatment in Mukhanov et al (1992) is less general than that of 
Bardeen (1980) and Brandenberger et al (1983) in two respects. First, they as- 
sume the anisotropic stresses are zero, and second, they make a speciflc choice of 
gauge invariants a priori, namely those associated with the so-called longitudinal 
gauge. 

Currently, increasingly accurate observations are driving theoretical cosmology 
towards more sophisticated models of matter and the study of possible nonlinear 
deviations from FL cosmology. Motivated by this state of affairs, our long term 
goal is to provide a general but concise description of nonlinear perturbations of FL 
cosmologies that will reveal the mathematical structure of the governing equations 
and enable one to make the transition between different gauge-invariant formula- 

^We follow the nomenclature of Wainwright and Ellis (1997) where an FL cosmology is a RW 
geometry that satisfies Einstein's field equations. 
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tions, thereby simplifying and relating the different approaches that have been used 
to datej^ In pursuing this objective we have found it necessary to revisit linear 
perturbation theory, even though it is by now a mature discipline!! Our intent in 
the present paper is to formulate the governing equations for the linear theory in a 
particularly simple and concise form in order to facilitate the extension to nonlinear 
perturbations. 

Based on earlier work by Bruni et al (1997) on gauge- invariant higher order per- 
turbation theory, Nakamura (2003) introduced a geometrical method for construct- 
ing gauge invariants for linear and nonlinear (second order) perturbations which he 
later applied to derive the governing equations (see Nakamura (2006) and Naka- 
mura (2007)). In the present paper we use a dimensionless version of Nakamura's 
method for constructing gauge invariants, but we complement it with the observa- 
tion that gauge invariants are of two distinct types: intrinsic gauge invariants, i.e., 
gauge invariants that can be constructed from a given tensor alone, and hybrid gauge 
invariants, i.e. gauge invariants that are constructed from more than one tensor. 

In Nakamura's approach, the linear perturbation of any tensor is written as the 
sum of a gauge-invariant quantity and a gauge-variant quantity, which is the Lie 
derivative of the zero order tensor with respect to a suitably chosen vector field X. 
A choice of X yields a set of gauge-invariant variables that are associated with a 
specific fully fixed gauge. We will show that for the metric tensor there exist two 
natural complementary choices of X that yield intrinsic metric gauge invariants. One 
choice, used in all of Nakamura's papers, leads to the two gauge-invariant metric 
potentials of Bardeen (1980), which are associated with the so-called Poisson gauged 
The other choice leads to the two gauge-invariant metric potentials of Kodama and 
Sasaki (1984), which are associated with the so-called uniform curvature gauged 
We will show that these two preferred choices lead to two distinct ways in which to 
present the linearized Einstein field equations: with the Bardeen choice the evolution 
of linear scalar perturbations is governed by a second order (in time) linear partial 
differential operator, while with the Kodama-Sasaki choice the evolution is governed 
by two coupled temporal first order linear operators. 

The plan of the paper is as follows. In Section [2] we discuss the geometrical 
construction of gauge-invariants: we focus on the metric tensor and, with the Ein- 
stein tensor and the stress-energy tensor in mind, on mixed rank two tensors. In 
Section [3] we use intrinsic gauge invariants to derive the general governing equations 
for linear perturbations in two gauge-invariant forms associated with the Poisson 
and the uniform curvature gauges. The required expressions for the Einstein gauge 
invariants are derived efficiently in Appendix [B], where we also give a general concise 
formula that expresses the Riemann gauge invariants in terms of the metric gauge 
invariants. One of the ingredients in our derivation is the so-called Replacement 
Principle, which is formulated in Appendix \M In Section H] we give an interpre- 

2 See, for example, Noh and Hwang (2004), Nakamura (2007) and Malik (2007). 

■^For some recent reviews and books, see, for example, Mukhanov (2005), Tsagas et al (2008), 
Weinberg (2008), Durrer (2008), Malik and Wands (2009) and Lytli and Liddle (2009). 

^The Poisson gauge, which was introduced by Bertschinger (1996) (see his equation (4.46)), is 
a generalization of the longitudinal gauge, which only applies to scalar perturbations. 

^ See, for example, Malik and Wands (2009), page 20, and other references given there. 
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tation of the intrinsic matter gauge invariants and specialize our equations to the 
cases of a perfect fluid and a scalar field. Section [5] contains a brief discussion of 
future developments. 



2 Geometrical definition of gauge invariants 
2.1 General formulation 

Following standard cosmological perturbation theory (see for example, Chapter 7.5 
in Wald (1984)), we consider a 1-parameter family of spacetimes gabi.^), where gab{0), 
the unperturbed metric, is a RW metric, and e is referred to as the perturbation 
parameter^ We assign physical dimension length to the scale factor a of the RW 
metric and (length)"^ to gab{^)- Then the conformal transformation 

gabie) = a^gabie), (1) 

yields a dimensionless metric gab{^)- Our reason for making this choic^ll concerning 
the allocation of physical dimensions is that it enables one to create dimensionless 
quantities by multiplying by the appropriate power of a, leading to simple pertur- 
bation equations that do not contain a explicitly. We refer to Appendix [Bl where 
this process is applied. 

The Riemann tensor associated with the metric gab{(^) is a function of e, denoted 
R"'^cd{^), as is the Einstein tensor, G^le). The stress-energy tensor of the matter 
distribution is also be assumed to be a function of e, denoted T%{e). We include all 
these possibilities by considering a 1-parameter family of tensor fields A[e), which 
we assume can be expanded in powers of e, i. e. as a Taylor series: 

A(e) = (0)A + eWA + ie2 (% + ... . (2) 

The coefficients are given bj|§ 



(°)A = A(0), «A , 

0€ 



(2)A 



(3) 



where ^^^A is called the unperturbed value, ^^^A is called the first order (linear) per- 
turbation and '•^^A is called the second order perturbation of A{e). 

The primary difficulty in cosmological perturbation theory is that the pertur- 
bations of a tensor field A{e) depend on the choice of gauge, and hence cannot be 
directly related to observations. It is therefore desirable to formulate the theory in 
terms of gauge-invariant quantities, i. e. to replace the gauge- variant perturbations 



^We use Latin letters a, &,...,/ to denote abstract spacetime indices. 

^An alternative choice in cosmology is to make a dimensionless and let the spacetime coordinates 
of 5ab(0) have dimension length (see, for example, Malik and Wands (2009), page 48). This choice 
is unsuitable for our purposes since it does not lead naturally to perturbative equations involving 
dimensionless quantities. 

^The notation A{e) should be viewed as shorthand for A{x, e), indicating that the tensor fields 
are functions of the spacetime coordinates, which necessitates the use of partial differentiation with 
respect to e. 
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(2)/^^ Qf yl(e) by gauge-invariant quantities. In this paper we restrict our 
attention to first order, i.e. hnear, perturbations, but with a view to subsequently 
working with higher order perturbations we use a method pioneered by Nakamura 
(2003), and adapt it so as to create quantities that are gauge- invariant and dimen- 
sionless. 

A linear gauge transformation is represented in coordinates by the equation 

5« = + + • • • , (4) 

where is an arbitrary dimensionless vector field on the background. Given a 
family of tensor fields A{e) the change induced in the first order perturbation ^^^A 
by a gauge transformation is determined by 

A(i)A = ^^("U, (5) 

where £^ denotes the Lie derivative with respect to and A^^^A := ^^^A — ^^^A (see, 
for example, Bruni et al (1997), equations (1.1) and (1.2)). We now introduce an as 
yet arbitrary dimensionless vector field X on the background which we use to define 
the dimensionless object^ 

(i)A[X] := a" ((^U - , (6) 

where we assume that A{e) is such that a^A{e) is dimensionless. It follows from ([5]) 
and ([6]) that 

A(i)A[X] = a"" (^^(OU - £ax^°^ = aJ'£^-^x^''^A. (7) 

The key step is to choose an X that satisfies 

AX'^ = C", (8) 

under a gauge transformation. With this choice, ([7]) implies that A^^)A[X] = 0, z.e., 
*^^^A[X] is gauge-invariant. We say that '^-'^)A[X] is the gauge invariant associated 
with ^^''A by X -compensation. Equations ([5]), ([6]) and ([8]) are central to our version 
of Nakamura's method for constructing gauge invariants associated with the first 
order perturbation of a tensor A (see Nakamura (2007), equations (2.19), (2.23) and 
(2.26)). In what follows we will drop the superscript '•^^ on A for convenience since 
in this paper we are dealing only with first order perturbations. 

The above 'gauge compensating vector field' X, which for brevity we shall refer 
to as the gauge field, requires comment. Unlike the geometric and matter tensor 
fields such as gabi^) and T%{e) it is not the perturbation of a corresponding quantity 
on the background spacetime. Instead it should be viewed as a vector field on 
the background spacetime that is constructed from the linear perturbations of the 
geometric and matter tensors in such a way that ([8]) holds. We will construct specific 
examples of X in section 12.21 We note that in choosing the gauge field X we are 
essentially fixing the gauge [i.e. making a choice of gauge), which is accomplished 

^When we consider second order perturbations in a subsequent paper, we will denote and 
X°' by and '■^•'X", and introduce a second pair of vector fields denoted '•'^^'^ and '-^^X". 
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in the traditional approach by making a choice of the vector field ^ that determines 
the gauge transformation^ One advantage of using the gauge field X is that one 
immediately obtains a geometric connection between the gauge invariants associated 
with different choices of gauge. This matter is discussed in more detail in Uggla and 
Wainwright (2011). 

Before continuing we briefly digress to point out that gauge invariants associated 
with a tensor A are of two distinct types: those that are solely constructed from 
components of ^^^A and ^^^A are called intrinsic gauge invariants, while those that 
depend on the components of another perturbed tensor are called hybrid gauge 
invariants. In particular if the gauge field X is formed solely from components of 
(^U and (°U, then A[X] is an intrinsic gauge invariant, otherwise A[X] is a hybrid 
gauge invariant. 

In the following sections we will calculate the quantities in equations ([5]) and ([6]) 
for various geometric objects A. To do this it is necessary to use the well known 
formulae for the Lie derivative. The formula for a tensor of type (1,1), which we 
now give, establishes the pattern: 

£^A% = A\,e + e,bA''c - C,cA%, (9) 

where ^ denotes partial differentiation. In a formula such as one can replace 
the partial derivatives by covariant derivatives. For our purposes it is convenient to 
use the covariant derivative °Va associated with the unperturbed conformal metric 

£^A% = (^,A%)^- + ("VfeDA'^e - (^cOA-'b. (10) 

We also need to work in a coordinate frame so that we can calculate time and 
spatial components separately. We thus introduce local coordinate^ = {rj,x^), 
with r] being the usual conformal time coordinate for the RW metric gabi^), and 
such that the unperturbed conformal metric jab '■= ^afe(O) has components 

7oo = -l, 7oi = 0, jij, (11) 

where jij is the metric of a spatial geometry of constant curvature. The curvature 
index of the RW metric, denoted K, determines the sign of the curvature of the 
spatial geometry, and if non-zero can be scaled to be +1 or —1 (see, for example, 
Plebanski and Krasinski (2006), page 261). 

The spacetime covariant derivative ^a determines a temporal derivative °Vo A = 
drjA, where 9^ denotes partial differentiation with respect to r], and a spatial covari- 
ant derivative that is associated with the spatial metric 7^. We introduce the 
notation 

:= '^iA. (12) 

^''See, for example, Malik and Wands (2009); equations (6.17), (7.3) and (7.4) provide an example 
in connection with the metric tensor. 

^^We use Greek letters to denote spacetime coordinate indices on the few occasions that they 
occur, and we use Latin letters i,j,k^m to denote spatial coordinate indices, which are lowered 
and raised using 7^ and its inverse 7*^ , respectively. 

^^Since we assigned a to have physical dimension length, the conformal time rj and the conformal 
spatial line-element ^ijdx^dx^ are dimensionless. We choose the to be dimensionless, which 
implies that the 7^ are also dimensionless. 
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The derivative operators dr, and Dj will be used throughout this paper once lo- 
cal coordinates have been introduced. However, for simplicity we shall denote the 
derivative of a function /(r^) that depends only on rj by f'{i]). 

With our present allocation of dimensions, the scalar H defined by 

n:=- = aH, (13) 
a 

where H is the Hubble scalar]^ is dimensionless. We shall refer to it as the dimen- 
sionless Hubble scalar. The use of this scalar, e.g. by Mukhanov et al (1992) (see 
page 218), is essential in eliminating a from the perturbation equations. 



2.2 Metric gauge invariants 

We expand gab{^), defined by equation ([1]), in powers of e: 

and label the unperturbed metric and (linear) metric perturbation according to 

7ab := ^'^-gab = -gabiO), fab ■= ^'^U = ^^(0), (14) 

which is consistent with (|3]). Applying the general transformation law ([5]) to the 
metric tensor gab{^) = o?gab{,^) we obtain 

A^^Vfe = £^''%ab, or, equivalently, Afab = a""^ £^{a^jab), (15) 

in terms of the notation f |T4l) . The gauge invariant fa6[^] associated with the metric 
perturbation fab by X- compensation, given by (EI), assumes the form 

fab[X] = fab - a~'^£x{a^lab)- (16) 

Introducing local coordinates and using ([9]) and ( |TOl) adapted to a (0, 2) tensor, 
equations (ITSl) and (fT6!) lead to 

A/oo = -2(a, + H)e, foo[X] = /oo + 2(a, + H)X\ (17a) 

A/o. = -D,e° + 5,6, fo.[X] = /oi + D,X° - a,X„ (17b) 

A/,, = 2H + 2D(,e,) , % [X] = f,, - 2HX°7,, - 2D(,X,) . (17c) 

In order to construct a gauge field X that satisfies ([8]), using only the metric, 
we need to decompose the metric perturbation fat into scalar, vector and tensor 

^•^Recall that H := where t is cosmic time, and that 4^ = a. 

a at ^ ' dr] 
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des0 We introduce the notatioio 



/oo = -2v^, (18a) 
/o. = D,E + S„ (18b) 
fij = + 2D,D,-C + 2D(,C,-) + 2^^, (18c) 

where the vectors Bi and Cj and the tensor Cij satisfy 

WBi = 0, WCi = 0, C\ = 0, U'dj = 0. 

The vector ^ is also decomposed into a scalar mode and a vector mode with com- 
ponents 

It follows from (dZD, ([IH]) and (dH) that 

A^={d, + n)e, AB = -e + dr,C, AC = C, Aij = -ne, (20a) 
AB, = a^e^, AC, = I,, (20b) 

ACij = 0. (20c) 

We can draw two immediate conclusions. First, it follows from fl20bl) and fl20cp 
that Bi — C[ and Cij are gauge invariants. We introduce the following bold-face 
notation: 

Bj := Bi — drfCi, Cij := Cij. (21) 
Second, by inspection of f|T9|) . fl20al) and fl20bl) we obtain 

A(D,C7 + Q) = e., Ax = A (^^^ = -^^ (22) 

where we have introduced the notation 



X:=B- d,C. (23) 

We are now in a position to satisfy the requirement (IE]). Firstly, referring to ( !22l) . 
we can satisfy the spatial part AX* = ^* of the requirement by choosing 

Xi = BiC + C\, (24) 

-^*In order to guarantee that the functions B,Bi,C,Ci and dj in (fT8|) are uniquely determined 
by /oi and /ij we need to assume that the inverses of D^, + 2K and + 3K exist. See the 
proposition in Appendix IB.ll See also Nakamura (2007), following equation (4.15), for a helpful 
discussion of this matter. 

^^We are denoting the scalar mode functions by ip, B, C and ip, in agreement with Mukhanov et 
al (1992) (see equation (2.10), but note the different signature) and Malik and Wands (2009) (see 
equations (2.7)-(2.12)), with the difference that we use C instead of E. Bardeen (1980) used the 
notation A, —B,Ht and — i?L + ^D^iJr for these functions, the choice of the fourth one being 
motivated by harmonic decomposition. Bardeen's notation has been used by subsequent authors, 
for example, Kodama and Sasaki (1984) and Durrer (1994), although the latter author replaced 
-B by B. 
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which we will take to be our default choice for Xi. With this choice, the ex- 
pressions ( !T7|) for the components of the gauge invariant fafe[X], when combined 
with ( ITSj) . assume the form 

foo[X] = -2$[X] , (25a) 

fo^X] =D,B[X]+B,,, (25b) 

fy[X] = -2^[X]7,, + 2Q,. (25c) 

where 

<^[X]:=^-{dr, + n)X^, ^X]:=^ + nX'^, B[X]:=x + X°, (25d) 



and Bj, Cjj and x ^^^^ given by fl2Tl) and fl23l) . respectively. 

Secondly, referring to f l22|) . we can satisfy the timelike part AX*' = of the 
requirement ()8]) in two obvious ways, by choosing 

X° = Xj:=-x, or X° = X° := (26) 

which leads to the metric gauge invariants associated with the Poisson gauge, or the 
uniform curvature gauge, respectively. On substituting these choices into ( I25d[) we 
obtain the conditions 

B[Xp] = and ^[Xe] = 0, (27) 
which characterize these two gauge choices. 

The Poisson gauge invariants 

On substituting the first of equations fl26l) into fl25l) we obtain 

foo[Xp] := -2$, foi[Xp]:=Bi, f^, [Xp] := -2^7^,- + 2Q,- , (28) 

where 

$ := <l>[Xp] = + (9, + 7^) X, ^ := ^[Xp] = ^-Ux. (29) 

Here $ and are the scalar metric gauge invariants associated with the Poisson 
gauge and ^ is the Bardeen potential. 

The uniform curvature gauge invariants 

On substituting the second of equations ( l26i) into ( l25i) we obtain 

foo[Xc] = -2A, fo^Xc] = DiB + Bi, f,,-[Xe] = 2Q„ (30) 

where 

A := $[Xc] = <^ + (a, + -H) ^, B := B[X,] = X - ^. (31) 

Here A and B are the scalar metric gauge invariants associated with the uniform 
curvature gaugejlll introduced by Kodama and Sasaki (1984)0 



^^The gauge-fixing conditions for the Poisson gauge are B ~ C — Q,Ci — Q va ([T! 
-^''The gauge-fixing conditions for tlie uniform curvature gauge are i/' = C* = 0, Ci = in p8| 
""^^See equations (3.4) and (3.5), noting that Hl + n~^HT = —"0 and B — k^^H!p = x- 
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In concluding this section we note that the gauge fields X used to construct the 
above gauge invariants have the same spatial components X* given by fl24|) in both 
cases, leading to ( 125|) . with the vector and tensor modes described by the gauge 
invariants Bj and Cij, respectively. The difference lies in the scalar metric gauge 
invariants which are related according tc0 

A = $ + (a, + ?^)^, B = -^, (32) 

as follows from fl29|) and fl3Tl) . In both cases the gauge invariants are intrinsic since 
the gauge field X depends only on the metric. 

A reader of this paper should be aware of the lack of agreement in the literature 
on labelling the scalar metric gauge invariants associated with the Poisson gauge. 
Our choice of in f l2I?]) is the one initiated by Mukhanov et al (1992), and 

subsequently used by Nakamura (see, for example, Nakamura (2006)) and Malik 
and Wands (2009). On the other hand Durrer (2008) and Liddle and Lyth (2000) 
reverse the roles and use (\I',$), while Kodama and Sasaki (1984) use (\I',— $). 
Bardeen's original notation is ($A) —^h)- 



2.3 Gauge invariants for mixed rank 2 tensors 

In this subsection we consider a rank two tensor A%, such that Aab is symmetric 
and a?A% is dimensionless. We expand A% in a Taylor series in e as in ([2]), and 
assume that '^^^A% obeys the background symmetries, which means it is spatially 
homogeneous and isotropic: 

D/°U'^;3 = 0, = (°U^o = 0, = I 5) (33) 

We introduce the notation 

A, := a^(-(°U°o + rA\), Cl := (34) 

where as before ' denotes differentiation with respect to r/. We further assume that 
A% satisfies the conservation law '^aA% = 0. It follows that in the background 

a\^'^A%)' = 3a'n{-^''>A'o + |WA\) = SHAa, (35) 

which, in conjunction with implies that 

A'a = -(1 + SCDHAa. (36) 

We can now calculate the gauge invariants A'fe[X] associated with ^^'^A% by X- 
compensation, as defined by equation with n = 2. It is convenient to decompose 
(^Mj- into its trace '■^^A'fc and tracefree part defined by 

(^U^- :=«A^--|«y4\5V (37) 

^^These relation have recently been given by Christopherson et al (2011). See their equations 
(4.22) and (4.23). 
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A straightforward calculation using (E]), ([9]), ( ITOl) and ( l33ll leads tc0 



A°o [X] = (iU°o - 371 AaX^ (38a) 

A° [X] = + ^^D,X°, (38b) 

A\ [X] = a" + mAAClX'' , (38c) 

A^-[X] = Wi^-. (38d) 



In deriving these equation we have used (IMl) and (1351) to express '•^^A'o, ^'^•'v4\ and 
their derivatives in terms of Aa and C^. 

Equation fl38dp implies that Aj[X] is an intrinsic gauge invariant since it is 
constructed solely from the components of ^^^A%. We denote this quantity by 

A):=A)[X] = a'm,. (39) 

One can form two additional intrinsic gauge invariants by taking suitable combina- 
tions of A°o[X], A°i[X] and A\[X]. Indeed it follows from ([38]) that 

A := ClA°o[X] + iA\,[X] = a\C'/''>A'o + 'i'^A%), (40a) 

A, := - {-DM[X] + SnA^lX]) = -a' (D,,(^U°o + 3^^^^°) , (40b) 

which implies that A and Aj are intrinsic gauge-invariants. 

In summary, the tensor A% can be described by the three intrinsic gauge invari- 
ants Aj, A, and Aj, given by f p9|) . f l40ap and fl40bp . and one hybrid gauge invariant 
A° [X], given by ( ]38b|) . In section [XT] we will use these objects, constructed in terms 
of the Einstein tensor and the stress-energy tensor, to give a concise derivation of 
the governing equations in gauge-invariant form for linear perturbations of PL. 



3 Linearized governing equations 
3.1 General formulation 

In this section we work with the linear perturbations of the Einstein tensor and the 
stress-energy tensor, denoted by ^^^G% and ^^^T%, and defined via equation (|3]). The 
corresponding unperturbed quantities are labelled by a superscript 

We begin by imposing the ba ckg round Einstein equations ^^^G% = ^^Xr%. The 
non-zero components are given bjo 

(o)G°o = -3(^2 + K) = - a2(o)p = a' (°>r°o, (41a) 

a2 Wg) = -(2H' + n^ + K)6) = a^^%6) = a^'-^yp), (41b) 

where is given by ( IT3l) and K is the curvature index. It follows from ( |4T1) . ( I34ll 
and ( l35l) . with A replaced by G and T, respectively, that 

Ag = 2{-H' + n^ + K), At = a\^% + (°)p), (42a) 

Aa = -{^ + ^Cl)nAG. ^T = ^r (42b) 

^°We do not include the ^^^A^q components since they can be expressed in terms of the other 
components and the metric perturbation, due to the assumed symmetry. 

^^See, for example, Mukhanov et al (1992), equation (4.2), noting the difference in signature. 
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The conservation law fl35p . with A replaced by T, gives 

a\(%y = -S-HAt = -ma\^% + (43) 

The background Einstein equations imply that Ac = At and Cq = C^. We denote 
the common values by A and C^: 

A = Ag = At, C^ = Cl = C|. (44) 

The linearized Einstein field equations are given by 

(1)^"^^ = ^^yr%. (45) 



In simplifying the linearized field equations we will make use of the intrinsic gauge 
invariants associated with the Einstein tensor and with the stress-energy tensor, 
which are given, in analogy with f l39|) . (140 ap and fl40bp . by 



= Wg) 




(46a) 


G. = -a' (D,(i)G'°o + 3H(^)G° ) , 


T, = -a' (D,(in + 3H(^>r° ) , 


(46b) 


G = a2(CS«G°o + |(^te\), 


T = a2(c|(i>rOo + |(^>r\), 


(46c) 


where 










(47) 


We also need the hybrid gauge invariants G*i [X] and T^i [X] , which are given 


by (|38b!) 


with A replaced by G and T: 






G°[X]=a2«GU^GD.X°, 




(48) 



Since the gauge invariants (H6|) and ( 148|) are linear in '■^-'G'^ and ^^\r% with co- 
efficients depending on ^^^G^ and '•''-T'^, respectively, it follows that the linearized 
Einstein field equations immediately imply the following relations: 

G) - f) = 0, G, - T, = 0, G - T = 0, (49a) 

G°[X]-T°[X] = 0. (49b) 

Expressions for the Einstein gauge invariants G^-, Gj, G and G^^-'^] in terms of the 
metric gauge invariants, decomposed into scalar, vector, and tensor modes, are given 
in equations flll6p and fll20p in Appendix [Bl To proceed we likewise decompose the 
matter gauge invariants TVTj,T and T'i[X] into scalar, vector, and tensor modes 
and label them as foUowso 



f) = B)U + 2f ^D(fcnj) + U), (50a) 

T, = DiA + A,, (50b) 

T = r, (50c) 

T°[X] =D,y[X] + T/„ (50d) 



In subsection 14. 1 1 we comment on the choice of the symbols 11, F, A and V. 
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where 

D^n, = 0, U\ = 0, D^n^' = 0, D^A^ = O, DV^ = O, (SOe) 

and 

:= D(,D,) - |7,,D2^ := D*D,. (50f) 

We stress that in making this decomposition we are not making any assumptions 
about the physical nature of the stress-energy tensor. By inspecting flll6p . fll20p 
and fl5U]) one concludes that equations (jl^ decompose into a scalar mode, a vector 
mode and a tensor mode, which we label as follows: 

BijA + D(iAj) + Aij = 0, 
DiB + B,; = 0, 
C = 0, 
DiE[X] + = 0. 

Since we are assuming that the inverses of the operators D^,D^ + 2K and + 
3K exist we can use the proposition in Appendix IB.ll to write the linearized field 
equations concisely as 

Scalar mode: A = 0, B = 0, C = 0, E[X] = 0. (51a) 
Vector mode: =0, B^ = 0, = 0. (51b) 

Tensor mode: Aij = 0. (51c) 



3.2 Scalar mode 

In this subsection we give the governing equations fl51al) for the scalar mode, first 
expressing them in terms of the uniform curvature gauge invariants A = $[A^c] and 
B = B[Xc] (see 0311) ). The scalars A, B and C in fl51ap are obtained without any 
calculation by taking the differences of equations (11161) and (ISUl) and reading off the 
scalar part. The scalar E,[X] is obtained in a similar manner from (I120p and (I50dp 
with X = Xp. The resulting equations ar^ 

(dr, + 2H) B + A = -n (52a) 

n [{dr, + BWjA + C^D^B] = |r + iD^n, (52b) 

-H + Sir) B = -iA, (52c) 

-HA + {\Ag - K)^ = -\V, (52d) 



where 



2%' 



^=^ + l + 3C^' (53) 

(see equation flll7p in Appendix |B]), and V = ^[-^p]- We shall refer to these 
equations as the uniform curvature form of the governing equations for the scalar 
mode. 
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In deriving (|52b|) we use (|52a| to replace {d^ + 2%) B + A by -11. 
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We now give the governing equations in terms for the Poisson gauge invariants \1/ 
and $. We ehminate A in (152bp using fl52al) and in (152d|) using (132!) . and ehminate 
B using = — The resulting equations are 

^ - $ = n, (54a) 

{c - c^D^) ^ = ir + (|d2 + nidr, + BH)) n, (54b) 

(D^ + 3K)^ = iA, (54c) 
dr^^ + -H^ = -\V, (54d) 

where the differential operator £ is defined by 

£(.) := -H^dr, + BH){d^ + 2H) (^) , (55) 

and B is given by ( 153|) . Expanding the brackets yieldH 

C = dl + 3{l + Cl) nd.^ + n^B - (1 + 3Cl)K. (56) 

We shall refer to the above equations as the Poisson form of the governing equations 
for the scalar mode, and to the evolution equation (]54bp as the Bardeen equation. 

Equations fl52|) and fl5^ . which are linked by the factorization property fl55l) . 
constitute one of the main results of this paper. Either system of equations determine 
the behaviour of linear scalar perturbations of an FL cosmology with arbitrary stress- 
energy content whose scalar mode is described by the gauge invariants F, 11, A and 
V . The structure of these two systems of equations differs in a significant way. 
In the system ( 152|) the time dependence is governed by two first order differential 
operators (9^ + BH and dr^ + 21-1, while in the system ( 15^ the time dependence is 
governed by the second order linear differential operator C A key point is that 
the coefficients in these operators depend only on the background KW geometry, and 
this dependence manifests itself through the appearance of 'H,'H','H" and K. This 
property is significant since it means that these operators will have the same form 
irrespective of the nature of the source in the FL background model, e.g. whether 
it is a perfect fluid with p = p{p), or a scalar field with potential V{(j)). What will 
differ, however, is the functional dependence of 7i{rj), which is determined by solving 
the Einstein equations in the background RW geometry, and hence depends on the 
source. Furthermore these differential operators will also appear in the linearized 
field equations in any geometrical theory of gravity, whose field equations depend in 
some way on the Einstein tensor. 

To the best of our knowledge equations (!52|) have not been given in the literature, 
although if one performs a harmonic decomposition one obtains a system of first 
order ordinary differential equations closely related to that given by Kodama and 
Sasaki (1984) (see Chapter 2, equations (4.6a-d)). Likewise, the governing equations 
in Poisson form f l5^ have not appeared in the literature in the above fully general 
form. The use of the Poisson gauge invariants was initiated by Bardeen (1980), and 



the evolution equation fl54bp for \E' is now commonly used, although it is written in 



a variety of different forms, as a partial or ordinary differential equation with the 



^Referring to (|42|) to express H' in terms of Ac and then use the equation for A'q. 
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coefficients usually expressed in terms of the matter variables of the background FL 
model. In contrast we have written the Bardeen equation in a fully general form 
in terms of the purely geometric differential operator C, which is defined by the 
factorization property ( 155|) . We can relate our form of the equation to the literature 
by expanding C as in (ISBj) and expressing the coefficients in terms of the matter 
variables. If the matter content is a barotropic perfect fluid and a cosmological 
constant and one imposes the background Einstein field equations then the geometric 
coefficients Cq and B can be written as 

Cl = cl n^B = {cl - w)pa^ + (1 + cl)Aa^ - (1 + 3cl)K, (57) 

using f l4T|) . fj44|) and f l79|) . The form in the literature that is closest to the purely 
geometric form fl56l) is that given by Mukhanov et al (1992), equation (5.22), who 
replace Cq by the matter quantity as in f lFTl) but retain 1-i and T-L' . Nakamura 
(2007) gives the same expression (see his equation (5.30)). A more common form in 
the literature has B, in addition to Cq, expressed in terms of the background matter 
variables as in ( 157|) . The earliest occurrence of which we are aware is Harrison 
(1967), equation (182), followed by Bardeen (1980), equation (5.30), after making 
the appropriate changes of notation and setting A = 0. See also Ellis, Hwang and 
Bruni (1989), equation (31) and Hwang and Vishniac (1990), equation (105) o 



3.3 Vector and tensor modes 



First, we give the governing equations fl51b|l for the vector mode. The vectors A, 



and Bj in fl51b|) are obtained without any calculation by taking the differences of 



equations (11161) and (!50|) and reading off the vector part. The vector Ej is obtained 
in a similar manner from (11201) and (]50dp . The resulting equations are 

{dr, + 2H)B, = -2Hi, (58a) 
(D^ + 2K)Bi = 2Vi, (58b) 



as well as the relation Aj = SHVi, which is satisfied identically (see equation 
If Hi is specified and can be regarded as a source term, the evolution equation fl58ap 
is a ffist order linear ordinary differential equation that determines Bj, which in 
turn determines Vi by differentiation using fl58bl) . 

Second, we give the governing equations fl51cp for the tensor mode. The ten- 
sor Ajj in f l51cp is obtained without any calculation by taking the differences of 
equations (I116P and (150|) and reading off the tensor part, leading to 

{d^^ + 2ndr, + 2K--D^)Cij = U,j. (59) 

If Hij is specified and can be regarded as source term, this is a second order linear 
partial differential equation that determines Cij. 



^^In these two references, the evolution equation in question arises in the 1 + 3 gauge-invariant 
approach to perturbations of FL, and the unknown is a vector quantity that is related to the scalar 
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4 Interpretations and examples 

4.1 Interpretation of the matter gauge invariants 

In this section we give the physical interpretation of the gauge invariants 11, F, A 
and V[X] associated with the scalar mode of the stress-energy tensor. 

We begin with the decomposition of a stress-energy tensor with respect to a unit 
timelike vector field u"", which is given by 

n = (p + P^Ub + p6% + {q'^Ub + u'^qb) + Ti\ (60) 

where 

M% = 0, ^"a = 0, Ua'K% = {]. (61) 

We choose to be the timelike eigenvector of T%, which implies = 0, i.e. we are 
using the so-called energy frame (see for example, Bruni et al (1992), page 37). 

Assuming that the unperturbed stress-energy tensor '^'^'>T% has the isotropy and 
homogeneity properties of the RW geometry, the expansion ([2]) to linear order for 
p,p, Ua and vr'fc has the form@ 

p = (o)p + e p = ^% + e (62a) 

= = 7r°i, TT^. = + e ^^Vj , (62b) 

uq = —a{l + e^p), Ui = a{0 + evi). (62c) 

Decomposing Vi into a scalar and vector mode yields 

V, = BiV + Vi, D'v, = 0. (63) 

We use boldface in writing Vj in view of the fact that this quantity is a dimensionless 
gauge invariant, as can be verified by applying ([5]) to -u^. 

For ease of comparison with other work, we note that the expansion of u"" = g°'^ui, 
to linear order, expressed in terms of f , v* and the linearly perturbed metric, is given 
by 

u"" = a~\l-e^), u' = a'^[Q + e{'D\v- B) + {y' - B'))]. (64) 

We digress briefly to mention that our expansion of the four- velocity differs from the 
usual approach in the literature in that we use the covariant vector Ua to define the 
perturbed three-velocity instead of the contravariant vector it°, since we find that 
this leads to a number of simplifications^ For example, Malik and Wands (2009) 
(see equation (4.4)) have 



u 



I ^[Q + e[T)'vMw + v'Mw)l 
so that 



VMW = V- B, Vljyy = V* - B\ 



^^The form of uq is determined by the requirement that is a unit vector. Recall that ip is one 
of the metric potentials in ^5)) . 

^^The source of these simplifications is the fact that Ui is invariant under purely spatial gauge 
transformations while is not. 
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From fl60l) and fl62p . and making use of ([3]), we obtain the following expressions 
for the components of the linear perturbation of the stress-energy tensor: 

= -(i)p, = 3 (^V, = {% + ^'^) = (^VV (65) 



It follows from (1461), (EQl) and ([65]), in conjunction with (I42j) and (143]), that the 
matter gauge invariants are determined by 



Ti)U + 2f '^D(fenj) + ff,-, (66a) 

r = a2(-C|«p+^'^p), (66b) 
A = a2((i)p+((o)p)'^), (66c) 

V[X] = At{v + X°), Vi = Ax^i. (66d) 

Before continuing we derive an additional relation. It follows from (I40b|) with A 
replaced by T that 

T, = -D,T°o [X] - mT\ \X\ . (67) 
On substituting from (I50bl) and (I50dp into this equation, we conclude that 

A = -T°o[X] - mV\X\, A, = -mVi. (68) 

We can now give the physical interpretation of the matter gauge invariants. 
First, the gauge invariants 11, Ilj and Iljj represent the anisotropic stresses. The 
interpretation of F is given in the context of a perfect fluid in the next section. 
Next, the gauge invariants V = V[Xp\ and Vi play a role in determining the shear 
and vorticity of Ua- The relevant formulae are given in (11281) in Appendix IB. 31 In 
particular, ^[Xp] determines the scalar mode of the shear according to 

BWj = |^r^D2(D2 + 3K)V[Xp], (69) 

as follows from fll28p in conjunction with ( I66dl) with X = Xp and the identity f ll26el) . 
We will hence use V := V"[Xp] as our standard choice for the gauge invariant 
However, since the choice V^[Xc] is also of interest we note that 

V[X,] - V[Xp\ = AtB, (70) 

as follows from f l66d]) . (126]) and (|3T]). 

Finally, in order to interpret A we need to make a small digression. For any scalar 
field A with the property that a'^A is dimensionless we can define a dimensionless 
gauge invariant A[X] according tc0 

A[X] = a" ((^U - ((°U )' X°) . (71) 

For the matter density p we denote the gauge invariant by p[X]: 

p[X] = a^(«p-(Wp)'X°). (72) 
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This is equation ([6]) specialized to the case of a scalar field. 
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On choosing X = Xy with X° := —v it follows from fl66cll that A = p[X^]. By 
comparing fl72|) with equation (3.13) in Bardeen (1980)o we conclude that p[X^], 
and hence A, equals the well-known Bardeen gauge-invariant density perturbation 
em, up to a factor of The specific relation is 

A = ia^%)em. (73) 

We note that the choice X^ = —v, in conjunction with our default choice (|24|) for 
the spatial components of X, is associated with the so-called total matter gauge 
(see, for example, Malik and Wands (2009), pages 23-24). Thus A is the density 
perturbation in the total matter gauge. In addition it turns out that A is closely 
related to the 1 + 3 gauge-invariant approach to perturbations of FL, pioneered by 
Ellis and collaborators (see for example, Ellis and Bruni (1989), Ellis et al (1989)), 
in which the spatial gradient of the matter density orthogonal to plays a key role. 
To elucidate the relation we define the dimensionless spatial density gradiento 

Vaie) = a^K\e) ^V,p(e), K\e) = 5^ + Ua{t)u\e). (74) 

A straight-forward calculation shows that 2^a(0) = and that to linear order 

(i)Do = 0, = D,A - mVi, (75) 

from which we conclude that A eq uals the the scalar mode of the linear perturbation 
of the spatial density gradient o In addition it follows from fl50b|) and fl68|) that 



Tj, giving a physical interpretation of the intrinsic gauge-invariant Tj 
To end this section we comment on our choice of notation. In using the symbols 
n, r, A and V for the matter gauge invariants we are following Kodama and Sasaki 
(1984) with the difference that we scale the variables as follows: 

Ii = a'pIiKs, T = a^pTKSi A = a^pAKs, V = AtVks, (76) 

where p and p refer to the background. Our choice of scalings simplify the equations 
considerably. 

4.2 Perfect fluid 

For a perfect fluid the matter gauge invariants are restricted according to 

n = 0, Hi = 0, U) = 0. (77) 

In addition it follows from ( I42bl) and (I66bl) that 

r = if and only if p = p{p), (78) 



^^One has to take into account differences in notation, the conservation equation P5)) . and the 
fact that Bardeen has performed a harmonic decomposition. 

■^°Our Va differs from that in Bruni, Dunsby and Eiiis (1992) by a factor of pa^ (see their 
equation (24)). 

3iNote that ' p{e) = /o(e)- 
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i. e. if and only if the equation of state is barotropic. In this case it is customary to 
introduce the notation 

cl := C|, w := (5^, (79) 

where = w if to is constant, as follows from (I42bp . 

On account of fl77j) the governing equations in the Poisson form fl5^ for scalar 
perturbations imply that \1' — $ = 0, which (in conjunction with the background 
field equations) reduces the governing equations for the scalar mode in the perfect 
fluid case to 

(£ - c^D^)^ = |r, (80a) 
(D^ + 3K)^ = |A, (80b) 
^' + = -\V, (80c) 

where C is given by ( 156|) with Cq = C^ = cj. and B is expressed in terms of the 
background matter variables according to (l57|l . 

4.3 Scalar field 

For a minimally coupled scalar field we show in Appendix [C] that the matter gauge 
invariants are given by 

r = (l-C|)A, (81a) 
V[X] = -^%'cj)[X], V, = 0, (81b) 

n = 0, Hi = 0, n^- = 0, (81c) 



where (j)[X] is the gauge invariant associated with (^^0 by X-replacement, given bjo 

= (82) 

Note that At and are given by ( llSSp . The governing equations (!54l) in Poisson 
form imply that \I' — $ = 0, and then reduce to 

(£-C'D2)^ = i(l-C2)A, (83a) 
(D2 + 3K)^ = 1A, (83b) 
^' + = i(°y(/)p, (83c) 

where (p^ := 0[Xp], and where we have used Cq = = C"^. By combining fl83ap 
and ( I83bl) we obtain an evolution equation for ^ without a source term: 

(£-3(1 -C2)K-D2) ^ = 0, (84) 

where £ is given by (1^ . Having solved this equation one can calculate 4)^ and A 
from fl83p . If one expresses in £ in terms of the unperturbed scalar field and its 
derivatives (see (11330 ) and sets K = 0, equation (|8^ coincides with equation (6.48) 
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This is a special case of equation (|7ip . 
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in Mukhanov et al (1992). For the generalization to arbitrary K, see Nakamura 
(2007), equation (5.39) S 

One can also use the governing equations fl52|) in uniform curvature form, ob- 
taining equations equivalent to those derived by Malik (2007) (see equations (2.20)- 
(2.23), noting that he is considering multiple scalar fields). 

5 Discussion 

We have given a systematic account of the gauge-invariant quantities that are as- 
sociated with a linearly perturbed RW geometry and stress-energy tensor, empha- 
sizing the role of intrinsic dimensionless gauge invariants. First, we have shown 
that there are two distinct choices of dimensionless intrinsic gauge invariants for the 
perturbed metric, which are the gauge invariants associated with the Poisson gauge 
and the uniform curvature gauge, through the work of Bardeen (1980) and Kodama 
and Sasaki (1984), respectively. Second, we have introduced dimensionless intrinsic 
gauge invariants for the Einstein tensor and the stress-energy tensor, which we used 
to derive a particularly simple and concise form of the governing equations for linear 
perturbations of FL models. The specific form of the governing equations for the 
scalar mode depends on the choice of intrinsic gauge invariants for the perturbed 
metric. The Kodama-Sasaki choice leads to a coupled system of two first order (in 
time) linear differential operators that govern the evolution of the uniform curvature 
metric gauge invariants (see equations ( 15^ ). On going over to the Poisson picture, 
the product of these two operators yields the second c»r(ier linear differential operator 
£ that governs the evolution of the Bardeen potential (see equation (!56|) ). thereby 
providing a link between the two forms of the governing equations. A common fea- 
ture of both systems is the appearance of the physically motivated gauge-invariant 
density perturbation A that is one of the intrinsic gauge invariants associated with 
the stress-energy tensor (see equations fl52cp and fl54cl) ). 

The mathematical structure of the governing equations for linear perturbations 
that we have elucidated here has in fact a much wider significance. Indeed, as 
one might expect on the basis of elementary perturbation theory, the governing 
equations for second order (nonlinear) perturbations have precisely the same form, 
apart from the inclusion of a source term that depends quadratically on the linear 
metric perturbation!^ As an illustration of this we give the form of the equations 
that govern second order scalar perturbations using the metric gauge invariants 
associated with the Poisson gauge: 



{c - c^D^) = i(2)r + (iD^ + nidr, + Bn)Y'^ + s,.oi{^'^), (85b) 



'We note a minor typo: a factor of 2 multiplying should be deleted. 

This behaviour has been noted in general terms by Nakamura (2006), equations (38)-(39). 







2^ ~l~ ^velocityi,^ '^) 
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where S',(^^¥) is a source term that depends quadratically on the first order gauge- 
invariant metric perturbation ^^'^ab = fab in equation ( 128!) . The key point is that, 
apart from the source terms, equations flSSp have the same form as equations fl54l) . 
with the variables ^'^^ and ^'^^^ being the metric gauge invariants at second order 
determined by the Nakamura procedure. The second order matter terms '^^•'T, 
^^^A and *^^V are defined in analogy with the first order terms 11, F, A and V 
after expanding the stress-energy tensor T% to second order in powers of e. All 
the complications lie in the source terms, whose explicit form has to be found by 
calculating the Riemann tensor to second order. In order to solve the above second 
order equations the source terms, which include scalar, vector and tensor modes, 
first have to be obtained by solving the governing equations for the scalar, vector 
and tensor linear perturbations. In a subsequent paper we will derive both the 
above Poisson form and the corresponding uniform curvature form of the governing 
equations for second order perturbations, relating our formulation to other recent 
work. 

In this paper we have focussed exclusively on using the linearized Einstein field 
equations to describe the dynamics of scalar perturbations. There are, however, two 
alternatives to the direct use of the linearized Einstein equations. First, one can 
use the linearized conservation equations for the stress-energy tensor, and second, 
one can use the 1 + 3 gauge-invariant formalism@ in which the evolution equations 
are obtained from the Ricci identities. An advantage of using the first approach 
independently of the Einstein equations is that the results are applicable to theories 
of gravity other than general relativity. An advantage of the second approach is that 
one initially derives exact nonlinear evolution equations, which are then subsequently 
linearized. Both of these approaches lead to a system of first order partial differ- 
ential equations that describe the evolution of scalar perturbations. An additional 
aspect of the dynamics of scalar perturbations that we have likewise not touched on 
in this paper is that under certain conditions {i.e. in the long wavelength regime) 
the governing equations admit so-called conserved quantities, i.e. quantities that 
remain approximately constant during a restricted epoch. These quantities, which 
are related to both the linearized Einstein equations and the linearized conservation 
equations, have been found to be useful in analyzing the dynamics of scalar pertur- 
bations during inflation. We refer to Uggla and Wainwright (2011), where we discuss 
the above aspects of the dynamics of scalar perturbations within the framework of 
the present paper. 
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A The Replacement Principle 

The expression for the perturbation of the Riemann tensor given in equation fll04p 
in Appendix [HI can be written symbohcally in the form: 

a^'-'^R^',, = I'^Uf), (86) 

where L" cd is a hnear operator and / is shorthand for fab- The Replacement Principle 
for the Riemann curvature states that the gauge invariants associated with 
and with fat by X-compensation are related by the same linear operator: 

K^UX] = l^UnX]), (87) 

where i[X] is shorthand for fa;,[X]. 

This result is adapted from more general results given by Nakamura (2005) (see 
in particular, his equations (3.12), (3.15) and (3.23)). Similar results hold for the 
Einstein and Weyl tensors. Use of the Replacement Principle in Appendix [B] makes 
the transition from gauge- variant to gauge-invariant equations particularly easy and 
transparent. 

B Derivation of the curvature formula 

In this appendix we derive expressions for the Einstein gauge invariants, namely, 
the three intrinsic gauge invariants Gj-,Gj and G, and the single hybrid gauge 
invariant '^^^G'ifX], defined by equations (146|1 and (HHjl . Our strategy incorporates 
the following ideas: 

i) Conformal structure. We adapt to the conformal structure of the background 
geometry, determined by the scale factor a of the RW metric, from the outset. 
In particular we create dimensionless quantities by multiplying with appropri- 
ate powers of a, which simplifies the equations considerably. 

ii) Index conventions. We represent tensors of even rank, apart from the metric 
tensor, with equal numbers of covariant and contravariant indices. This makes 
contractions trivial to perform and ensures that the components of the tensor 
have the same physical dimension as the associated contracted scalar. 

iii) Timing of specialization. We defer performing the decomposition into scalar, 
vector and tensor modes as long as possible, and do not make harmonic de- 
compositions. This strategy helps to reveal structure in the equations and 
serves to reduce the amount of calculation. 

Calculation of i?"^cd(e) 

We begin by deriving an exact expression for the Riemann tensoiQ R"'\d{^) of the 
metric Qabi^) in terms of the covariant derivative of the conformal background metric 



'We use the sign convention of Wald (1984) for defining the Riemann tensor. 
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■Jab- We thus relate the covariant derivative of gab{^), denoted ''Va, to that of jab = 
(7afe(0), denoted °Va. The relation is given by an object = Q^b defined by 

Q%c = g'^'Qdbc = Ig"" {^cgdb - "^dgbc + ^^bgcd) , (88) 

(see Wald (1984) equation (D.l)), with the property thatlf^ 

'VaA\ = "^aA^ + Q\dA% - Q'^acAV (89) 

It is convenient to write Q%i. as the sum of two parts: 

Q\c{t) = Q\c{t) + Q\c{e). (90) 

First, the transformation from ^Vq to "^Vq, which is associated with the conformal 
transformation gabi,^) = o?gab{,^), is described by 

Q%,{t) = 25%r,) - r\e)-g,,{e)u, (91) 

wherj^ 

:=OV,(lna) (92) 

(see Wald (1984), equation (D.3)). It follows that ^a^b = ^bra- Second, the trans- 
formation from'^Va to °Va, the covariant derivatives associated with (jabi^) and gab{0), 
respectively, is described by 

Q%c{e) = ir'ie) {"Vcgdbie) - ^.c(e) + °V, ^,,(e)) . (93) 

It follows from ^ajbc = that 

QV(0) = 0. (94) 

To calculate R"'^cd{^) we first perform the conformal transformation from gab to 
(jab, which yields 

a^R'^'M = + ^6\u\{e), (95) 

where 

U\{e) = - [i' (^V, - r,) + \5\r^rf] r^, (96) 

and R'^^cdi^) is the curvature tensor of the metric ^ab(e) (see Wald (1984), equation 
(D.7)). Second, by performing the transition from'^Va to °Va we obtain 

R^U^) = g''R\cd{e) = {^R\,d + 2°V[,Q"d]e + 2g'^/[cQ^d]e) , (97) 

where '^R°'bcd is the curvature tensor of the metric jab (see Wald (1984), equation 
(D.7)). The term 2g^^%M]e in can be written afl 

2^'^ %CQ\ = (°V[C Q\fW .gg^ 



■^^This example establishes the pattern for a general tensor. 

■^*Note that we always use the vector Tq in covariant form, since Va is independent of e, whereas 



r°- = g'^''{e)ri, is not. 
^^^Note that ^R^Kd = l''" "R^ecd 
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which we use to rearrange f l97|) . in conjunction with the relation ^cQ""^ 
In summary, R"'\d{^) is given by equation fl95l) with 



[c V|e| gd]f - Te/S- 



df 



re, 



where we have used "'VaVb = ^a^b — Q^abfc in obtaining f l99bp from fl96|) . 



-2Q('^''l 

(99a) 
(99b) 



Calculation of ^^^R" 



cd 



We now calculate the perturbation '^^^R"'\d of the Riemann tensor, defined via equa- 
tion (E]), expressing it in terms of the covariant derivative °Va associated with 
and the metric perturbation fab = ^^^Qab (see ([ED). We note that 



{l)gab 



(100) 



where the indices on are raised using_7"''. It follows from (|3]), f l93|) f l95|) and fl99|) . 
in conjunction with (EH) and ffTOOD . thaO 



a 



2(1)^. 



'''cd = ^''^R''\d + ^5^''\Pu\ 



where 



b 
abc 



I (°^c /afe — °Va fbc + °Vfe /ca) • 



bd 



Introducing local coordinates = (rjjX^) as in section [2?T] leads to 



(lOla) 

(lOlb) 
(lOlc) 
(lOld) 

(102) 



In addition we note that the quantity ^R"'bcd, the curvature tensor of the metric jab, 
is zero if one index is temporal, while if all indices are spatial 



^R% 



km 



(103) 



where the constant K describes the curvature of the maximally symmetric three- 
space. Equation (1101 p . in conjunction with (I102p and (I103p . yields the following 
expressions: 



wherj^ 



a2«i?%^ = |[D^D„ + in' - 'H')6Ufoo + (5, + n)Y^n,, (104a) 

a2Wi?%„ = 2D[,y^^], (104b) 

a'^'^R%m = -2 (D[fcD[' + K6i^') + 4?^(J[jY^fl, (104c) 

y,, = '^%,nfoo - D(,/,)o + ia,/.,. (io4d) 



40 



Note that R'^^cdie) depends on e through gab{(),g°'''{e) and Q'^ab[el)- 
iNote that Q% = -B^Jj^o + \flj- 
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Calculation of the Riemann gauge invariants 

We now apply the Replacement Principle to fll04p . which entails performing the 
following replacements: 

U-^fab[X], Y,,~^Y,,[X], a^^'^R'^id-^K'^UX], (105) 

where the gauge invariants are defined by equation (l6l). All components of the 
Riemann tensor can be obtained from the 'curvature spanning set' {R%j, R^^jk, 
R^^jm) or, alternatively, their spatial traces and their trace-free parts: 

/ jDOm pOm pfcm \ / D^i £>0i r)im \ ^'1nfi^ 

{-TL Om, -rt jm, -rt km), {-n- Oj; jk, ^ jm), [LVD) 

where 

oOi oOi 1 ri pOm him pirn 1 ri pfcm {Tr\^f^\ 

n Qj — n Oj — ^OjU Om, ^ jm — ^ jm ~ ^"jU km, l^iU/aj 

= — S\kR^"^j]m- (107b) 

Our motivation for choosing these particular components as the spanning set is that 
the first set of terms in fll06p are invariant under spatial gauge transformations, 
while the hatted quantities are fully gauge-invariant, as follows from 

We denote the gauge invariants associated with the spanning set f ll06p by 

l-Tt Omi^J, ■K' jm[-^\, ■> km[-^\), Oj , -K- jk, ^ jm), UUSJ 

and refer to them as the Riemann gauge invariants. As indicated by the notation 
{i.e. no dependence on the gauge field X) the hatted quantities are intrinsic gauge 
invariants. We now substitute the expression^ for fab[X] given by ( |25l) into the 
bold- face version of (11041) . and calculate the gauge invariants (11081) . It is convenient 
to split Yij into a trace and a trace- free part: 

Y,, = Y,, - l7,,Y, Y = Y\, (109) 

and to use the trace- free second derivative operator Djj defined in (150fp . We obtaicQ 



R'^^omix] = - [D^ + 3{n' - n')] ^x] + (a, + n) y[x], (iioa) 

R%,- = -D^$[X] + {dr, + n) Y)[X], (110b) 

R'^'kmlX] = 4 [(D^ + 3K) ^X] + nY[X]] , (110c) 

TV'^j^ = B)^[X] + -HY^iX] - (D^ - 2K) C), (llOd) 

R^'^.^IX] = |D,Y[X] - D„Y™ [X], (llOe) 

R°^-fc = 2D[,Y\,][X] + D^Y"^y[X]5\], (llOf) 



^^In using these expressions we are making the choice for Xi given in equation (|24p . Choosing 
Xi in this way simphfies the calculation but not the final form of the Riemann gauge invariants, 
since, as mentioned earlier, the spanning set is invariant under spatial gauge transformations. 

''^Use the identities (p^Bc)) . (fngd)) and (fT26h)) . 
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where 

Y[X] = -3{dr,^X] + n^X]) - T>^B[X], (llOg) 
Y,j[X] = -D,,-B[X] - D(,B,) + drid,. (llOh) 

These equations constitute one of the main results of this paper. They express the 
Riemann gauge invariants (11081) in terms of the metric gauge invariants fl2^ . They 
depend only on the choice of the temporal gauge field as can be seen from (I25d[) . 



Calculation of the Einstein gauge invariants 

The Einstein tensor and the Weyl conformal curvature tensor are defined in terms 
of the Riemann tensor according to 

G% := R% - \ b\R, where R\ := R := R\, (111a) 

C"^'cd := R''\d - 2 5\ + I ^\ ^\ R. (111b) 

The curvature spanning set (]106p can be replaced with the following spatially irre- 
ducible components of the Einstein tensor and the Weyl tensorQ 

['^ 0, m, '^i, {'^ Oj) jk), l-LJ-^J 

where 

• ^3 3 J 



:= G) - U^G^^. (113) 



It follows from f lTTT]) that 

— ~2-^ ^'^^ ^ ra — ~2\^ km + Omj, (^ii4aj 

(-^ j — -fL fm, Lr j — K oj + K jm, [114:0) 

O Oj — -^[K Oj — n jra), ^ jk — -tl jk- [LL4C) 

The Einstein gauge invariants, as defined by equations fl3^ . fl40ap and fl40bp with 
A replaced by G, can be expressed in terms of the curvature spanning set fll06p by 
using the bold-face version of f lll4p . This yields 

G\- := G^^] = ^% + (115a) 
- (D,G°o[X] + 3HG° [X]) = |D,R'='"fc„[X] - 3?^R°'"„„[X], (115b) 
C'aG%[X] + lG^m[X] = -I ((1 + 3C2 )R^™fc„[X] + 4Rn^[X]) . (115c) 

We find that it is simplest to express the Einstein gauge invariants (11150 in 
terms of the uniform curvature metric gauge invariants A and B defined by ( |3T1) . 
We accomplish this directly by choosing X = Xc in flllOp . and noting that by ( l27p 
we have ^I/[Xc] = 0. After simplifying using the identities f ll26ep and (I126fp we 
obtairS 



'j 
G, 

G 



Gij = BijG - D(i {d^ + 2n) Bj) + {d^ + + 2K- D^) dj, (116a) 
Gi = 2n'Di{'D^ + 3ir)B + |7/(D2 + 2K)Bi, (116b) 
G = 2'H[id^ + BWjA + 4 D'B] - fD^G, (116c) 



''^Note that C*4m = -'^C%^k i^^'L] 

in an orthonormal frame 



**^Here for convenience we use = 7ifeG 
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where we have introduced the notation 

OO/' 

G := -[A + {dr, + 2n)B], B:= — + l + SC^ . (117) 
We also need 

G°[X]=R0-^[X]. (118) 

We choose X = Xp in this equation, and using f lllOp in conjunction with the iden- 
tity (]126f[) we obtain 

G° [Xp] = -2D, (9^^ + 7^$) + I (D^ + 2K) B,. (119) 

We now use ( l32i) to express the right side of this equation in terms of A and B, 
which yields 

G%[Xp] = -2Di {HA + {^Ag - K)B) + ^ (D^ + 2K) B^. (120) 
The Weyl tensor 

The perturbation of the Weyl tensor is automatically gauge-invariant on account of 
the Stewart- Walker lemma (Stewart and Walker (1974)) since the Weyl tensor is zero 
in the background. We thus use bold-face notation for its components. From f llMcp 
we obtain 

Oj — a U oj — gl-Tt Oi ~ -K. jmj, ^ jk — 0, U jfe — ±t jk- U^J-J 

The Weyl tensor has a simpler form if we use Poisson gauge invariants and hence 
we choose X = Xp in (11 101) . Noting that B[Xp] = leads to 

C%, = -i [D^(^ + $) + d^B) - {d^ + - 2K) C)] , (122a) 
C°jfc = — 2D[j (B\] — (9^C\]) — JDmB'^[j6\i, Bij := DjjBj]. (122b) 

B.l Uniqueness of the decomposition into modes 

Proposition: If the inverses of the operators D^, + 2K and + 3X exist, then 
the equation 

Bi = BiB + S„ with B'ISi = 0, (123) 
determines B and Bi uniquely in terms of Bi, and the equation 

= D,,C + D(,C,) + A„ (124) 

with 

D*Cj = 0, Cij = Cji, C\ = 0, D*Cjj = 0, 

determines C, Cj and Cij uniquely in terms of Cij. In particular, ii Bi = then 
B = 0, Bi = 0, and if dj = then C = 0, d = 0, Cij = 0. 

Proof. Apply D* to f ll23p obtaining D*i?j = D^B. Using the inverse operator of 
this equation determines B, and then f ll23p determines Bi uniquely in terms of Bi. 
Next, apply D*-'' and D* to (HMD, obtaining 

B'^Cij = |d2(d2 + 3K)C, B'C^j = |Dj(D2 + 3K)C + (D^ + 2K)Cj. (125) 

By using the inverse operators these equations, in conjunction with fll24p . succes- 
sively determine C, Cj and Cij uniquely in terms of Cij . □ 
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B.2 Identities 

In obtaining our results we found the following identities useful: 





KS\iAj], 


(126a) 


D[fcD^]A^^' = 


2K6[k^'A^p, 


(126b) 


4(D[,D[^ + U = 


(D,^ + 1 (D^ + 3K) 40 ^, 


(126c) 


4 (D[fcD[^ + = 


(D2 - 2K)C\, 


(126d) 


BjU^A = 


|D,(D2 + 3K)A 


(126e) 


D^D(,A,) = 


i(D2 + 2K)A,, 


(126f) 


BiB^A' = 


(D2 + 2K)DiA\ 


(126g) 


Si^'A r'l = 




(126h) 



where Aij = Aji, Cij = Cji, C\ = and DjCj = 0. 
B.3 Kinematic quantities 

The kinematic quantities associated with a timelike congruence m° are defined by 
the following decomposition into irreducible parts: 

VaUb = -Uailb + H{gab + Ma^fe) + C^afe + ^^ab- (127) 

A routine calculation starting with equations (162|) - (!M|) and (!89|) applied to Ua yields 
the following non-zero components: 

a^^^H = [lT)\v - x) - (<9^^ + -Hip)] , (128a) 

Ui := = Bi{ip+ {dr, + V)v) + (9^ + WjWi, (128b) 

a) := a^%) = W,{v - x) + f 'D(fc (v,) - B,)) + 9,C^-, (128c) 

a;^-:=a«a;^-=Y'=D[fcV,], (128d) 

with the bold-face quantities being gauge-invariant on account of the Stewart- Walker 
lemma. 

C Scalar field 

A minimally coupled scalar field (j) is described by a stress-energy tensor of the form 

T% = VyVb(f> - [I V>V,0 + U{cP)] n, (129) 

with the associated Klein-Gordon equation V^Vc^ — t/^ = 0, where the potential 
U {(p) has to be specified. This stress-energy tensor is of the form fl60l) with 

p + p= -VcPVaCp, p-p = 2U{(P), 7Cab = 0. (130) 

When evaluated on the RW background, equation fll30p leads to 

a2((o)p + (o)p) = ((oy)2^ (0)^ _ (o)p ^ 2U{%). (131) 
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On using fll3ip to calculate ^'^^p', the conservation equation (H3!l leads to 

b" + 2n%' + a^U^ = 0, (132) 



which is the Klein- Gordon equation in the RW background. Further, by means 
of (112]), (03]), ffml) and (fT32|) we obtain 

1.1^.4(1^1^). (133) 

Viewing T% and (f) as functions of the perturbation parameter e, we can use (11290 . 
in conjunction with ([3]), to calculate ^^^%, obtaining 

(i>f^- = 0, (^>rO = -(oy D,(iV, (i>rOo + i(^>r^, = -2[/,/^)0. (134) 



It follows using ( l38ll with A replaced by T and (I133p . that the matter gauge invariants 
assume the form 

f) = 0, T",[X] = -(°VD,(/>[X], n[X] + l^aX] = ~2a%^cf>[X], (135) 
where is the gauge invariant associated with by X-replacement, given by 

= (^V- (136) 

Equations (I135p and ( 150|) immediately lead to the expressions for the matter gauge 
invariants ( ISlbp and (I81cp . including 



= (137) 

Equation (1135^ . in conjunction with (1133^ and (1137^ . yields 

T%[X] + lT\[X] = -3(1 - Cl.)nV[X]. (138) 

We now substitute (11381) into the expression for F given by (I46cp and (I50cp to 
obtairQ 

F = (1 - C|)(-T°o[X] - 3nV[X]). (139) 
which on comparison with fl68l) leads to equation (18 lap . 
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